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Abstract 

We provide analytical solutions to the thermodynamic Bethe ansatz 
equations in the large and small density approximations. We extend 
results previously obtained for leading order behaviour of the scaling 
function of affine Toda field theories related to simply laced Lie alge- 
bras to the non-simply laced case. The comparison with semi-classical 
methods shows perfect agreement for the simply laced case. We derive 
the Y-systems for affine Toda field theories with real coupling constant 
and employ them to improve the large density approximations. We 
test the quality of our analysis explicitly for the Sinh-Gordon model 
and the {G^ , D± )-affine Toda field theory. 
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1 Introduction 



In the context of 1+1-dimensional integrable quantum field theories numer- 
ous methods have been developed to compute various quantities in an exact 
manner, that is non-perturbative in the coupling constant. Sometimes it is 
even possible to perform the related computations analytically. Often the 
evaluation within one particular approach lacks information, typically a con- 
stant, which might be supplied by an entirely different method. Ideally one 
would like to achieve a situation in which each approach is self-consistent. 

An example for the situation just outlined is given for instance in the form- 
factor program []]], which allows in principle to compute correlation functions. 
The lowest non-vanishing form-factor is not fixed within this approach and 
is typically obtained from elsewhere. For instance the vacuum expectation 
value of the energy-momentum tensor can be extracted from the thermody- 
namic Bethe ansatz (TBA) ||. Alternatively, one can compute correlation 
functions by perturbing around the conformal field theory ||. Also in this 
approach one appeals to the thermodynamic Bethe ansatz for the vacuum 
expectation value of the energy-momentum tensor and to the Bethe ansatz 
0] for the relation between the coupling constant and the masses. The latter 
correspondence is also needed in an approach initiated recently in |§, 01, 
where it was observed that the ultraviolet asymptotic behaviour for many 
theories may be well approximated by zero-mode dynamics. The consider- 
ations in || exploit the knowledge of an exact reflection amplitude, which 
on one hand results from certain manipulation on the three-point function 
of the underlying ultraviolet conformal field theory and on the other hand 
has a semi-classical counterpart in the related quantum mechanical problem. 
The question why this approach allows to compute scaling functions with a 
high accuracy is still to be settled [§]. 

The computation of scaling functions by means of the TBA does con- 
ceptually not require any additional input from other methods. However, 
up-to-now it is only possible to tackle the problem numerically due to the 
nonlinear nature of the central equation involved. Several attempts have 
been made to formulate analytical approximations. This is desirable for var- 
ious reasons, one being that the numerical effort becomes quite considerable 
for some models with increasing particle species content. A further reason is 
of course that analytical expressions allow to study further the deeper struc- 
tures of the theories. For instance for affine Toda field theories (ATFT) [[J 
related to simply laced Lie algebras some analytical expressions have been 
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provided []T0] , [H, [12] . In the approximation method of [|TI| [11], 0] a constant 
was left undetermined, which can be fixed in the same spirit as outlined for 
the other methods in the preceding paragraph, namely by appealing to an- 
other approach. Contrary to the claim in @, we demonstrate in the present 
manuscript that it is possible to fix the constant in this way without ap- 

also to 
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proximating higher order terms. Extending the analysis of [|Tl 
the non-simp ly laced case in the present manuscript, we will demonstrate 
in addition that the constant may also be well approximated from within 
the TBA-analysis. Furthermore, we give simple analytical expressions for 
improved approximations in the large and small density regime. 

Our manuscript is organized as follows: In section 2 we provide large 
and small density approximations for the solutions of the TBA-equations. In 
section 3 we assemble the necessary data for ATFT needed to extend our 
previous analysis to the non-simp ly laced case. We derive universal TBA- 
equations and Y-systems for all ATFT and show how they may be utilized 
to improve on the analytical approximations. We derive the related scaling 
function. We test the quality of the various approximations for the explicit 
example of the Sinh-Gordon model and the (G 2 
conclusions in section 4. 



(i) n( 3 ) 



DD-ATFT. We state our 



2 Large and small density approximations 
2.1 The TBA 

The object of investigation of the TBA is a multiparticle system containing 
n different particle species with masses rrii (1 < i < n) whose dynamical 
interaction is described by a factorizable diagonal scattering matrix Sij(9), 
which is a function of the rapidity difference 9. We assume the statisti- 
cal interaction to be of fermionic typeQ. Adopting the notation of [I2| the 



thermodynamic Bethe Ansatz equations 0, which characterize the thermo- 
dynamic equilibrium of such a system, are the n coupled nonlinear integral 
equations 

n 

rrrn coshfl + In (e il(0) - l) = ^ * Lj) (9) . (1) 

i=i 

*In order to keep the discussion as simple as possible we do not treat general statistics 
here as for instance Haldane type | fl3[ . Generalizations of our arguments in this sense are 
straightforward. 
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Here the scaling parameter r is given by the inverse temperature times a 
mass scale. The convolution of two functions is abbreviated as usual by 
(/ * g )(p) ■= l/2?r / d9'f(6 - B')g{9'). The TBA kernel reads 

^(d)--=-^S^6). (2) 

The functions Lj, which are to be determined as solutions of the TBA equa- 
tions ([!]), are related to the particle densities p\ and the densities of available 
states p l h as Lj = ln(l+p*/pjj, such that for physical reasons Lj > 0. Keeping 
this definition in mind, we speak of the large density regime when Lj > In 2 
and of the small density regime when Li < In 2. It is sometimes useful to 
express matters in terms the pseudo-energies := — ln[exp(Lj(0)) — 1)]. 

Having solved the TBA-equations ([!]) for the L-functions one is in prin- 
ciple in the position to evaluate the scaling function 

oo 

c(r) = ^ J2 m * / ddL i( 6 ) cosh6 ( 3 ) 
i=i { 

which can be interpreted as off-critical effective central charge belonging to 
the conformal field theory obtained in the ultraviolet limit, i.e. r — * 0. It is 
our goal in this manuscript to approximate this function in a simple analytical 
way to high accuracy. 



2.2 Approximative analytical Solutions 

In general it is possible to solve the TBA-equations numerically, where the 
convergence of the iterative procedure is guaranteed by means of the Banach 
fixed point theorem However, the numerical problem becomes quite 



complex when one increases the number of particle species. For this reason, 
and more important because one would like to gain a deeper structural in- 
sight into the solutions of ([!]), it is desirable to obtain analytical solutions to 
the TBA-equations. Due to the nonlinear nature of (p]) only few analytical 
solutions are known. Nonetheless, one may obtain approximated analytical 
solutions when r tends to zero. For large Li, i.e. for large particle densities, it 
was shown in [[10], 11, that the integral equation ([!]) may be turned into a 



set of differential equations of infinite order. Under certain natural assump- 
tions, which are however not satisfied universally for all models, one may 
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approximate these equations by second order differential equations, whose 
solutions are given by 



L°{0) = In ( W CW for \9\ < aroC06 ^. (4) 



The restriction on the range of the rapidity stems from the physical require- 
ment Li > 0. The n constants r] { = ^ . 77^ are determined by a power series 
expansion of the TBA kernel 

/oo °° 
de ^(e)e^ = 2^H)^t n ■ (5) 
■°° n=0 

The dependence on the scaling parameter r enters through the quantity 

7T 

A = 2(tf 4 -ln(r/2)) ' (6) 

Here the 5, are constants of integration. There is a very crude lower 
bound we can put immediately on <5j. From the fact that L®(0) > In 2, we 
deduce <5j > l/7r/ v /^~ + ln(r/2). For particular models we will provide below 
a rigorous argument which establishes that in fact they do not depend on 
the particle type, such that we may replace f3 i — > j3 and Si — > 5. We will 
also show that they can be fixed by appealing to the semi-classical approach 
in || [/[. In addition, we provide an argument which determines them ap- 
proximately from within the TBA analysis by matching the large and small 
density regimes. Since the constant turns out to be model dependent, we 
will report on it in detail below when we discuss concrete theories. 

The restriction on the range for the rapidities in for which the large 
density approximation ceases to be valid, makes it desirable to develop 

also an approximation for small densities. For extremely small densities we 
naturally expect that the solution will tend to the one for a free theory. 
Solving ([I]) for vanishing kernel yields the well-known solution 

L{{9) = ln(l + e- rm!COshe ) . (7) 

Ideally we would like to have expressions for both regions which match at 
some distinct rapidity value, say to be specified below. Since L®{6) and 
L{(9) become relatively poor approximations in the transition region between 
large and small densities, we seek for improved analytical expressions. This 
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is easily achieved by expanding ([I]) around the "zero order" small density 
approximations. In this case we obtain the integral representation 



^(^) = exp(-rm l coshe + ^" =i (^.*Lj)(^)) . (8) 

For vanishing ip^ we may check for consistency and observe that the functions 
L s {9) become the first term of the expansion in (0). One could try to proceed 
similarly for the large density regime and develop around L\ instead of L{. 
However, there is an immediate problem resulting from the restriction on 
the range of rapidities for the validity of L°, which makes it problematic to 
compute the convolution. We shall therefore proceed in a different manner 
for the large density regime and employ Y-systems for this purpose. 

In many cases the TBA-equations may be expressed equivalently as a 
set of functional relations referred to as F-systems in the literature [14[ . 
Introducing the quantities Yj = exp (—£,), the determining equations can 
always be cast into the general form 

Y0 + ^^{9 - i-Kfi) = exp(0i(0)) (9) 

with n being some real number and gi(9) being a function whose precise 
form depends on the particular model. We can formally solve the equation 
by Fourier transformations 

Y0) = exp [( 9i * 7i )(0)] , 7i (0) = [2/icosh(0/2///)] 1 / 2 (10) 

i.e. substituting (|ToD into the l.h.s. of (||) yields exp (g^O)). Of course this 
identification is not completely compelling and we could have chosen also a 
different combination of Vs. However, in order to be able to evaluate the 
gi(6) we require a concrete functional input for the function Yi(6) in form of 
an approximated function. Choosing here the large density approximation 
L° makes the choice for gi(6) with hindsight somewhat canonical, since other 
combinations lead generally to non-physical answers. 

We replace now inside the defining relation of gi(6) the Y's by Yi (8) — > 
exp(L°(6 1 )) — 1. Analogously to the approximating approach in JTD|, [TT|, [T2|], 
we can replace the convolution by an infinite series of differentials 

em = -(9i*im =-5> (m) ^w , (ii) 

771=0 
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where the z/s are defined by the power series expansion 



/oo oo _ 

The £7 m denote the Euler numbers, which enter through the expansion 1 / cosh x 
X]™ =0 x 2m i?2 m /(2m)!. In accordance with the assumptions of our previous 
approximations for the solutions of the TBA-equations in the large density 
approximation, we can neglect all higher order derivatives of the L®(9). Thus 
we only keep the zeroth order in (pp . From fll2"|) we read off the coefficient 
= 1/2, such that we obtain a simply expression for an improved large 
density approximation 

L\(9) = ln[l + Yl(6)] = ln[l + exp(fc(0)/2)] . (13) 

In principle we could proceed similarly for the small density approximation 
and replace now Yj (8) — > exp(Lf (#)) — 1 in the defining relations for the g^s. 
However, in this situation we can not neglect the higher order derivatives of 
the L\ such that we have to keep the convolution in ([□]) and end up with 
an integral representation instead. We now wish to match L\ and L\ in the 
transition region between the small and large density approximations at some 
distinct value of the rapidity, say Q™ . We select this point to be the value 
when the function fi(6) = (6/7r 2 )rmjZ/j(#) cosh#, which is proportional to 
the free energy density for a particular particle species, has its maximum in 
the small density approximation 



. (14) 



In regard to the quantity we wish to compute, the scaling function (|3]), this 
is the point in which we would like to have the highest degree of agreement 
between the exact and approximated solution, since this will optimize the 
outcome for c(r). Having specified the 0™, the matching condition provides 
a simple rational to fix the constant Si 

L\m = mr) =► c • (15) 

Clearly, in general we can not solve these equations analytically, but it is 
a trivial numerical problem which is by no means comparable with the one of 
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solving ([I]). Needless to say that the outcome of ( |T5| ) is not to be considered 
as exact, but as our examples below demonstrate it will lead to rather good 
approximations. One of the reasons why this procedure is successful is that 
L|(#™) is still ver Y close to the precise solution, despite the fact that is at its 
worst in comparison with the remaining rapidity range. 

Combining the improved large and small density approximation we have 
the following approximated analytical L-functions for the entire range of the 
rapidity 
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,L\{9) !<>.' 9 < (>■ 

Ld6) = {Lm for \9\>9? ' (16) 



such that the scaling function becomes well approximated by 

oo 

n „ 

c(r) * J dB ^) ■ 



(17) 



To develop matters further and report on the quality of L°, U , L s , L l we 
have to specify a particular theory at this point. 



3 Affine Toda field theory 

Affine Toda field theories form a well studied class of relativistic inte- 
grable quantum field theories in 1 + 1 space-time dimensions. To each of 
these field theories a pair of affine Lie algebras (Xn\x^) |15| is associated 
whose structure allows universal statements concerning its properties, like 
the S-matrix, the mass spectrum, the fusing rules, etc. Here X^ denotes a 
twisted affine Lie algebra w.r.t. a Dynkin diagram automorphism of order £. 
Both algebras are chosen to be dual to each other, i.e. X® is obtained from 
the non-twisted algebra Xn of rank n by exchanging roots and co-roots. For 
Xn simply-laced both algebras coincide, i.e. X n — X^ £ \ £ — 1, which is 
reflected in the quantum theory by a strong-weak self-duality in the coupling 
constant. Moreover, the mass spectrum renormalises by an overall factor 
and the poles of the S-matrix in the physical sheet do not depend on the 
coupling constant. For non-simp ly laced Lie algebras these features cease to 
be valid. The quantum masses are now coupling dependent and flow between 
the classical masses associated with Xn and X^ in the weak and strong 
coupling limit, respectively. Consequently, the physical poles of the S-matrix 
shift depending on the coupling and the strong-weak self-duality is broken. 
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3.1 The universal S-matrix 



Remarkably, despite these structural differences the S-matrix of ATFT can 
be cast into a universal form covering the simply-laced as well as the non- 



simply laced case [0, |I7 |. For our purposes the formulation in form of an 



integral representation is most useful 

S l3 {9) = exp / ^.(t)sinh-, (18) 
Jo r m 

hfi) = 8sinh(M,)sinh(t,^t) ([K] q{t)m )^ . (19) 

Denoting a q-deformed integer n as common by [n] q = (q n — q~ n )/{q l — q~ r ), 
we introduced here a "doubly q-deformed" version of the Cartan matrix K 
T6|, [I8|, K7\ of the non-twisted Lie algebra 



[KaU = (qq U + q^q'VSij ~ [hjh (20) 

for the generic deformation parameters q, q. The incidence matrix 1^ = 
25ij — Kij of the xjp related Dynkin diagram is symmetrized by the integers 
U, i.e. Iijtj = Ijiti. With ojj being a simple root we fix the length of the 
long roots to be 2 and choose the convention ti = Inside the integral 

representation fll8|) we take 



q{t)=e td \ q(t)=e MH , with h := , V H := (21) 

for the deformation parameters, where < B < 2 is the effective coupling 
constant. We further need the Coxeter numbers h, h and the dual Coxeter 
numbers h w , h w of Xn and X^\ respectively, as well as the £-th Coxeter 
number H = £h of X^. Complete tables of these quantities for individual 
algebras may be found in |Tjj| . 



The incidence matrix satisfies the relation 16, 17 



5^ [ Jy]?(«r) m 3 - = 2 cosh(0 fc + t^) , (22) 
i=i 

which will turn out to be crucial for the arguments below. We introduced 
here the imaginary angles 9h = infih and 9h = itt'&h- 
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3.2 The TBA-kernel 



From the universal integral representation fllq), we can now immediately 
derive the Fourier transformed TBA-kernel @ for ATFT. However, when 
taking the logarithmic derivative one has to be careful about interchanging 
the derivative with the integral, since these two operations do not commute. 
Comparison with the block representation of the S-matrix W% yields 



1 f°° t9 
( PM = -7T dt^exp-, (23) 

such that the Fourier transformed universal TBA-kernel (^) acquires the form 
<Pij(t) = -K^ij^t) = -87T ^DhtTT&hSiDhtjTr&H {[K]q(rt)q(irt)) ti ■ (24) 



To be able to carry out the discussion of the previous section we require the 

(2) 

1 y 



second order coefficient r)\j in the power series expansion (^). From (21) we 
read off directly 



Vlf = 7^7 B{2 - B)KJ tj = B(2 - B) (A, • A,) . (25) 



In the latter equality we used the fact that the inverse of the Cartan matrix 
is related to the fundamental weights as \ = J2j ^ij 1 a j> U = and 
H = £h = £h v . This implies on the other hand that 

Vt = ^B(2-B)(X l .p) (26) 
with p = At being the Weyl vector. Therefore 



1=1 



We used here the Freudenthal-de Vries strange formula p 2 = h v /12 dim X^ 1 ' 
(see e.g. ||19|| ) and the fact that dimXn 1 '' = n(h+l). Thus we have generalized 
the result of |12[ to the non-simply laced case. Notice that in terms of 
quantities belonging to the non-twisted Lie algebra Xn the formula ( p7| ) is 
identical for the simply laced and the non-simply laced case. 
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3.3 Universal TBA equations and Y-systems 



In analogy to the discussion for simply-laced Lie algebras ||12|| , the universal 
expression for the kernel ( 24l) can be exploited in order to derive universal 
TBA-equations for all ATFT, which may be expressed equivalently as a set 
of functional relations referred to as F-systems. Fourier transforming (H) in 
a suitable manner and invoking the convolution theorem we can manipulate 
the TBA equations by using the expression (|24l) . After Fourier transforming 
back we obtain 

n n 
3=1 3=1 

The universal TBA kernels A and T are then given by 

lM = ( 2 (^ + t ^H)cosh w J^y 1 , (29) 

r«(0) = ^iie + ^k-i-i^es), (30) 

k=l 

= [li{0+{0 h -t i B H )) + li {e-{d h -t i B H ))]8 ij . (31) 

The key point here is that the entire mass dependence, which enters through 
the on-shell energies mjcosh^, has dropped out completely from the equa- 
tions due to the identity (|22|). Noting further that 



[hj\q{i-w) m j cosh# = rrij cosh [9 + (2k — 1 — I^Oh] , (32) 

k=l 

we have assembled all ingredients to derive functional relations for the quan- 
tities Yi = exp(— Ei). For this purpose we may either shift the TBA equations 
appropriately in the complex rapidity plane or use again Fourier transforma- 
tions, see |12| 

Yl (e+e h + ue H ) y % (e-e h - ue H ) = [1+ n y '^ +eh ' t ' ejf)][1+r ' (e ' eh+ ^ Jf)] • (33) 

ft ff [i+Y- 1 (e+(2k-i~h 3 )e H )] 

3=1 fc=l 

These equations are of the general form @ and specify concretely the quanti- 
ties \x and gi{6). We recover various particular cases from (^). In case the as- 
sociated Lie algebra is simply-laced, we have Oh + UOn — ► iir/h, Qh — ^fiu — * 
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m/h(l — B) and — > 0, 1, such that we recover the relations derived in [12[ . 
As stated therein we obtain the system for minimal ATFT [[14) by taking the 
limit B — > zoo. 

The concrete formula for the approximated solution of the Y-systems in 
the large density regime, as defined in (|T3|), reads 



1 ii i 



yi(p\ - cos(2e/3 i )+cos(2(e h -f I e g )/3 i ) rrrr f 2^/3f \a 



j=l m=l 



Exploiting possible periodicities of the functional equations ( p3[) they may 
be utilized in the process of obtaining approximated analytical solutions f2(| . 



As we demonstrated they can also be employed to improve on approximated 
analytical solution in the large density regime. In the following subsection we 
supply a further application and use them to put constraints on the constant 
of integration 5i in (||). 

3.4 The constants of integration (3 and 5 

There are various constraints we can put on the constants P t and Si on general 
grounds, e.g. the lower bound already mentioned. Having the numerical data 
at hand we can use them to approximate the constant. In |12|] this was done 



by matching L° with the numerical data at 9 = and a simple analytical 
approximation was provided <5 num = ln[5(2 — B)2 1+B ^ B ^]. Of course the 
idea is to become entirely independent of the numerical analysis. For this 
reason the argument which led to ([T5|) was given. 

When we consider a concrete theory like ATFT, we can exploit its partic- 
ular structure and put additional constraints on the constants from general 
properties. For instance, when we restrict ourselves to the simply laced case it 
is obvious to demand that the constants respect also the strong-weak duality, 
i.e. Pi(B) = Pi(2 - B) and 81(B) = 5^2 - B). 

Finally we present a brief argument which establishes that the constants 
Pi are in fact independent of the particle type i. We replace for this purpose 
in the functional relations fl33|) the Y- functions by Y^{9)) and consider the 
equation at 9 = 0, such that 



cosh 2 [7r/^ fe + UDh)} - 2p 2 lVl _ A t4 
cosh 2 [vrA(^ + ^)] fiBt 



(cosh 2 [nPj (2m-l-I ZJ )■& g] \ 
cosh[7r/3 i (2m-l-/i j )i? H ] 

(35) 
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Keeping in mind that j3 i is a very small quantity in the ultraviolet regime, 
we expand ([35|) up to second order in /3 i , which yields after cancellation 



= d B ^~ B ) = ^K.RHi = B(2-B) T Jl 

j—l I I j—l I I 

(36) 

We substituted here the expression (E6l) for the constants t\a in the last equal- 



ity. Using once more the relation A, = Y2j a j > we can evaluate the inner 
product such that fl31f) reduces to 



n n nl 
j=l k=l Pi 

Clearly this equation is satisfied if all the /? f are identical. From the unique- 
ness of the solution of the TBA-equations follows then immediately that we 
can always take {i i — > (3. Since the uniqueness is only rigorously established 
12|j for some of the cases we are treating here, it is reassuring that we can 



obtain the same result also directly from (0). From the fact that the f5 i are 
real numbers and all entries of the inverse Cartan matrix are positive follows 
that 0\ = 0* for all i and j. The ambiguity in the sign is irrelevant for the 
use in L{6). 

3.5 The Scaling Functions 



In [E| it was proven that the leading order behaviour of the scaling function 



is given by 



377 / 7T 2 fi(2 - B)(h + 1) . 

c ( r ) " n ~ rx i / /o^^ ? =n[l- -rh ] \ „, • (3? 



(<5-ln(r/2)) 2 V. 4/i(<J - ln(r/2)) 2 

From our arguments in section 3.2, which led to the general expression for 
the constant 77 in form of (^7|), follows that in fact this expression holds for all 
afline Toda field theories related to a dual pair of simple affine Lie algebras 
(Xn\x^). However, strong- weak duality is only guaranteed for £ = 1. 

Restricting ourselves to the simply laced case, we can view the results 
of [|, |^, 0] obtained by means of a semi-classical treatment for the scaling 
function as complementary to the one obtained from the TBA-analysis and 
compare directly with the expression (B3). Translating the quantities in 
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|5|, [7j to our conventions, i.e. R — > r, B — > 5/2, we observe that c(r) 
becomes a power series expansion in (3. We also observe that the second 
order coefficients precisely coincide in their general form. Comparing the 
expressions, we may read off directly 

H *r(i-i)r(i + S)J-'" (39) 

for all ATFT related to simply laced Lie algebras[|. Here 7 £ denotes Euler's 
constant and k = (Yl^i^)^ 1 * s a constant which can be computed from 
the Kac labels n; of the related Lie algebra. Contrary to the statement made 
in M, this identification can be carried out effortlessly without the need of 
higher order terms. Recalling the simple analytical expression # num of [0 
we may now compare. Figure 1 demonstrates impressively that this working 
hypothesis shows exactly the same qualitative behaviour as 5 semi and also 
quantitatively the difference is remarkably small. 

To illustrate the quality of our approximate solutions to the TBA-equations, 
we shall now work out some explicit examples. 

3.6 Explicit Examples 

To exhibit whether there are any qualitative differences between the simply 
laced and non-simply laced case we consider the first examples of these series. 

3.6.1 The Sinh-Gordon Model 

The Sinh-Gordon model is the easiest example in the simply laced series and 
therefore ideally suited as testing ground. The Coxeter number is h = 2 in 
this case. An efficient way to approximate the L-functions to a very high 
accuracy is 



In 



-, cos(2/36»)+cosh(7r/3(l-B)) 
1 I ZT^3 



for 101 < 9 



ni 



L a (9) = { L ' ^ J 1 1 " (40) 

\exp [-rm cosh 9 + (cp * Lf){6)) for |0| > 9^ 



^The expressions in || and only coincide if in the former case m = 1 and in the 
latter m = 1/2. In addition, we note a missing bracket in equation (6.20) of ||, which is 
needed for the identification. Replace C — > — 4QC therein. 
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with 




' 0,0 0,5 1,0 Q 1,5 2,0 



Figure 1: Numerically fitted constant <5 num versus the constant from the semi- 
classical approach for the Sinh-Gordon model <J semi . 

The determining equation for the matching point reads 

sinh# m - rm/2sinh(2# m ) + cosh(fl m )(</ * L f )(9 m ) = . (42) 

For instance for B = 0.4 this equation yields 8 m = 11.9999 such that the 
matching condition ([T^) gives 5 m = 0.4913. Figure 2(a) shows that the 
large and small density approximation L° and U may be improved in a 
fairly easy way. In view of the simplicity of the expression L a the agreement 
with the numerical solution is quite remarkable. Figure 2(a) also illustrates 
that when using the constant <5 semi instead of 5 m the agreement with the 
numerical solutions appears slightly better for small rapidities. When we 
employ # num instead of 5 semi the difference between the two approximated 
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solutions is beyond resolution. However, as may be deduced from Figure 
2(b), with regard to the computation of the scaling function the difference 
between using 5 m instead of <5 semi is almost negligible. Whereas in the former 
case the resulting value for the scaling function is slightly below the correct 
value, it is slightly above by almost the same amount in the latter case. More 



on the approximation of the scaling function in form of ( |38"D may be found 
in 
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3.6.2 (Gcp, D^)-ATFT 

In this case we have h = 6 and H = 12 for the related Coxeter numbers. The 
two masses are mi = m sin(7r(l/6 — B/24)) and m 2 = m sin(7r(l/3 — 5/12)). 
The L-functions are well approximated by 



Lt(6) 



Mi + ^^g^^ v^^y 1 for ^ *™ 

exp[— rm\ cosh^ + ((p n * L[ + (p 12 * L{){6)} for \6\ > Q\ 



ln[1 cos(2^) + cos h ( :/ 3(l-||)) ^ / y 1 f |e| ^ 

L 4r; 2 /3 2 V cos 2 (/3(0+-|§)) J 11-2 



exp[— rm 2 cosh^ + (cp 21 * L{ + (p 22 * L 2 )(9)} for |#| > 6 2 



mi,/ ■ lr2] ^ 1 - - Tf\fa\\ \a\ ^ a™ 

with ip given by (^3|) and 

5^5(2-5) _ 7^(2 - B) _ 7**B(2 - B) 

m- 72 , V2- g , v- 36 • l«; 

Using now the numerical data L\(0) = 4.2524 and /^(O) = 3.67144 as bench- 
marks, we compute by matching them with L1(0) and L 2 (0) the constant 
to 5 = 1.1397 in both cases. This confirms our general result of section 3.4. 
Evaluating the equations (|15|) and (|T4] ) we obtain for B = 0.5 the matching 



values for the rapidities 9? = 12.744 and 6% = 12.278 such that 5? = 1.9539 
and #2™ = 1.5572. Figure 2(c) and 2(d) show a good agreement with the 
numerical outcome. 

The approximated analytical expression for the scaling function reads 

C(r) " 2 - 12(^-ln(r/2)r ^ 



This expressions differs from the one quoted in [12|], since in there the sign 
of some scattering matrices at zero rapidity was chosen differently. 
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4 Conclusions 

We have demonstrated that it is possible to find simple analytical solutions 
to the TBA-equation in the large and small density regime, which approx- 
imate the exact solution to high accuracy. By matching the two solutions 
at the point in which the particle density and the density of available states 
coincide, it is possible to fix the constant of integration, which originated in 



the approximation scheme of pi], [11], [L2| and was left undetermined therein. 
Alternatively the constant may be fixed by a direct comparison with a semi- 
classical treatment of the problem. It is not necessary for this to proceed to 
higher order differential equations as was claimed in |5[]. Of course one may 
proceed further to higher orders, but since the solutions to the higher order 
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differential equations may only be obtained approximately one does not gain 
any further structural insight and moreover one has lost the virtue of the 
first order approximation, its simplicity. 

We derived the Y-systems for all ATFT and besides demonstrating how 
they can be utilized to improve on the large density approximations we also 
showed how they can be used to put constraints on the constant of integra- 
tion. 

We have proven that the expression fl38|) for the scaling function is of a 
general nature, i.e. valid for all ATFT. It is desirable to extend the semi- 
classical analysis also to the non-simply laced case. This would allow to 
read off the constant 5 also in that case. 

Acknowledgments: The authors are grateful to the Deutsche Forschungs- 
gemeinschaft (Sfb288) for financial support. We acknowledge constructive 
conversations with B.J. Schulz at the early stage of this work and are grate- 
ful to V.A. Fateev for kind comments. 



References 

[1] P. Weisz, Phys. Lett. B67 (1977) 179; 

M. Karowski and P. Weisz, Nucl. Phys. B139 (1978) 445. 

[2] ALB. Zamolodchikov, Nucl. Phys. B342 (1990) 695; 

T.R. Klassen and E. Melzer, Nucl. Phys. B338 (1990) 485; B350 (1991) 
635. 

[3] ALB. Zamolodchikov, Nucl. Phys. B348 (1991) 619. 

[4] V.A. Fateev, Phys. Lett. B324 (1994) 45; 

ALB. Zamolodchikov, Int. J. Mod. Phys. A10 (1995) 1125. 

[5] A.B. Zamolodchikov and ALB. Zamolodchikov, Nucl. Phys. B477 (1996) 
577. 

[6] C. Ami, C. Kim and C. Rim, Nucl. Phys. B556 (1999) 505. 



17 



[7] C. Ahn, V.A. Fateev, C. Kim, C. Rim and B. Yang, Reflection Amplitudes of 
ADE Toda Theories and Thermodynamic Bethe Ansatz, |hep-th / 9907072) . 

[8] ALB. Zamolodchikov, 4-th Bologna Workshop on CFT and Integrable Models, 
(Bologna, 1999). 

[9] A.V. Mikhailov, M.A. Olshanetsky and A.M. Perelomov, Commun. Math. 
Phys. 79 (1981) 473; 

G. Wilson, Ergod. Th. Dyn. Syst. 1 (1981) 361; 

D.I. Olive and N. Turok, Nucl. Phys. B257 [FS14] (1985) 277. 

[10] ALB. Zamolodchikov, "Resonance Factorized Scattering and Roaming Tra- 
jectories", Preprint ENS-LPS-335 (1991). 

[11] M.J. Martins, Nucl. Phys. B394 (1993) 339. 

[12] A. Fring, C. Korff and B.J. Schulz, Nucl. Phys. B549 (1999) 579. 

[13] A. Bytsko and A. Fring, Nucl. Phys. B532 (1998) 588. 

[14] Al.B. Zamolodchikov, Phys. Lett. B253 (1991) 391. 

[15] V.G. Kac, "Infinite Dimensional Lie Algebras" (CUP, Cambridge, 1990). 

[16] T. Oota, Nucl. Phys. B504 (1997) 738. 

[17] A. Fring, C. Korff and B.J. Schulz, Nucl. Phys. B567 (2000) 409. 

[18] E. Frenkel and N. Reshetikhin, Deformations of W-algebras associated to 
simple Lie algebras, q-alg/9708006 v2. 

[19] P. Goddard and D.I. Olive, Int. J. Mod. Phys. Al (1986) 303. 

[20] Al.B. Zamolodchikov, Nucl. Phys. B358 (1991) 497; Nucl. Phys. B358 
(1991) 524; Nucl. Phys. B366 (1991) 122. 



18 



